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We  describe  a  Cs  atomic  magnetometer  operating  in  the  spin-exchange-relaxation-free  (SERF)  regime.  With 
a  vapor  cell  temperature  of  103  °Cwe  achieve  intrinsic  magnetic  resonance  widths  A B=  17  p G  corresponding 
to  an  electron  spin-relaxation  rate  of  300  s-1  when  the  spin-exchange  rate  is  Lse=  14  000  s-1.  We  also  observe 
an  interesting  narrowing  effect  due  to  diffusion.  Signal-to-noise  measurements  yield  a  sensitivity  of  about 
400  pG/ yHz.  Based  on  photon  shot  noise,  we  project  a  sensitivity  of  40  pG/yHz.  A  theoretical  optimization 
of  the  magnetometer  indicates  sensitivities  on  the  order  of  2  pG /  yHz  should  be  achievable  in  a  1  cm3  volume. 
Because  Cs  has  a  higher  saturated  vapor  pressure  than  other  alkali  metals,  SERF  magnetometers  using  Cs 
atoms  are  particularly  attractive  in  applications  requiring  lower  temperatures. 
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I.  INTRODUCTION 

Sensitive  atomic  magnetometers  have  recently  found  ap¬ 
plication  in  the  field  of  magnetic  resonance  imaging  [1], 
magnetoencephalography  [2],  and  searches  for  physics  be¬ 
yond  the  standard  model  [3].  A  recent  review  may  be  found 
in  Ref.  [4].  The  most  sensitive  atomic  magnetometers  pres¬ 
ently  are  the  spin-exchange-relaxation-free  (SERF)  magneto¬ 
meters  [5]  in  which  relaxation  due  to  spin-exchange  colli¬ 
sions  is  eliminated  by  operating  in  the  regime  where  the 
spin-exchange  rate  is  much  greater  than  the  rate  of  Larmor 
precession  [6,7].  In  Ref.  [5],  sensitivity  of  5  pG/  \/Hz  was 
achieved  with  the  cell  operating  at  190  °C  using  potassium 
atoms.  Estimates  of  the  fundamental  sensitivity  limit  of  such 
magnetometers  are  several  orders  of  magnitude  better  for  a 
1  cm3  volume  and  scale  as  the  square  root  of  the  spin- 
destruction  cross  section. 

Here  we  demonstrate  the  operation  of  a  Cs  magnetometer 
in  the  SERF  regime,  achieving  a  sensitivity  of  about 
400  pG/yHz  with  a  vapor  cell  temperature  of  only  103  °C 
The  overlapping  volume  of  the  pump  and  probe  beams  is 
about  0.02  cm3,  but  the  effective  volume,  determined  by  dif¬ 
fusion,  is  about  1  cm3.  Based  on  optical  rotation  measure¬ 
ments,  the  projected  photon  shot  noise  limit  for  our  experi¬ 
mental  conditions  is  about  40  pG/  y  Hz.  The  spin-destruction 
cross  section  for  Cs-Cs  collisions  is  ctsd=2X10-16  cm2, 
about  200  times  larger  than  for  K  [8].  Hence,  the  fundamen¬ 
tal  sensitivity  of  a  Cs  SERF  magnetometer  should  be  roughly 
a  factor  of  14  worse  than  that  of  a  K  SERF  magnetometer. 
However,  it  is  often  the  case  that  environmental  noise  due  to, 
for  example,  Johnson  currents  in  magnetic  shields  is  far 
larger  than  atomic  shot  noise,  and  hence  little  is  lost  by  using 
Cs.  One  of  the  primary  motivations  for  investigating  Cs  in 
the  SERF  regime  is  that  Cs  has  the  highest  saturated  vapor 
pressure  of  all  the  stable  alkali  metals,  yielding  significantly 
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lower  operating  temperatures.  This  opens  up  the  possibility 
of  operating  in  the  SERF  regime  with  low-buffer-gas- 
pressure  paraffin-coated  cells,  of  interest  because  this  allows 
operation  with  significantly  reduced  light  power.  The  lower 
temperature  required  for  a  Cs  magnetometer  is  also  attractive 
for  applications  such  as  NMR  measurements  with  liquids  in 
microfluidic  channels,  which  is  expected  to  be  an  important 
measurement  modality  in  future  “lab-on-a-chip”  devices  [9]. 
Cs  is  also  of  interest  for  a  search  for  a  permanent  electric 
dipole  moment  because  it  is  a  heavy  atom  where  CP  violat¬ 
ing  effects  are  enhanced  [10]. 

II.  BLOCH  EQUATIONS 

A  full  treatment  of  the  system  requires  the  use  of  density 
matrix  theory  (see,  for  example.  Ref.  [11]).  However,  the 
description  can  be  greatly  simplified  when  the  spin-exchange 
rate 

rsE=^E  =  «^SE^  (1) 

(here  n  is  the  alkali-metal  number  density,  crSE~2 
X  10-14  cm-2  is  the  spin-exchange  cross  section,  and  v  is  the 
average  relative  velocity  of  the  colliding  alkali-metal  atoms) 
is  much  faster  than  precession  in  the  magnetic  field, 
Fse^> gsfjLBB / (21+  l) .  Here  gs~ 2  is  the  electron  Lande  fac¬ 
tor,  pB  is  the  Bohr  magneton,  B  is  the  magnitude  of  an  ap¬ 
plied  magnetic  field,  and  I  is  the  nuclear  spin.  In  this  case  the 
density  matrix  assumes  a  spin-temperature  distribution  and 
the  ground  state  can  be  well  described  by  Bloch  equations 
for  the  electron  spin  polarization  P  =  (S)/.S’  [12-14]  as  fol¬ 
lows: 

f  =  °V2P  +  ^P)  X  »  +  R(s  -  P)  -  -  rpiP)  ' 

(2) 

Here  D  is  the  diffusion  coefficient,  s  is  the  optical  pumping 
vector  along  the  direction  of  propagation  of  the  pump  with 
magnitude  equal  to  the  degree  of  circular  polarization,  R  is 
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FIG.  1.  (Color  online)  Experimental  setup.  A  four-layer  mag¬ 
netic  shield  is  employed.  Circularly  polarized  light  tuned  to  the  D 1 
line,  propagating  in  the  z  direction  produces  ground  state  orienta¬ 
tion  in  the  z  direction.  The  x  component  of  orientation  Sx  is  detected 
via  optical  rotation  of  linearly  polarized  light,  tuned  to  the  wing  of 
the  pressure  broadened  D2  line,  propagating  in  the  x  direction. 


the  optical  pumping  rate  due  to  the  pump  beam,  and  rpr  is 
the  rate  of  depolarization  due  to  the  linearly  polarized  probe 
beam.  In  Eq.  (2),  Tl  and  7’2  are  the  relaxation  times  appro¬ 
priate  for  components  of  the  polarization  parallel  and  trans¬ 
verse  to  B,  respectively.  The  quantity  q(P )  is  the  nuclear 
slo wing-down  factor,  which  for  nuclear  spin  7=7/2,  is  [15] 


q(P)  = 


2(T,6+  17T,4  +  35T,2+  11) 
P6  +  IP4  +  lP2+\ 


(3) 


In  the  low  polarization  limit  q(0)  =  22,  while  in  the  high  po¬ 
larization  limit  </(l)  =  8.  The  latter  limit,  when  all  atoms  are 
pumped  into  the  stretched  state,  corresponds  to  the  slowing- 
down  factor  for  nuclear  spin  7=7/2  in  the  absence  of  spin- 
exchange  collisions  <7 = 27  +1  =  8.  The  transverse  relaxation 
time  can  be  written 


1 

T~i 


Tsd  + 


(4) 


where  rSD  is  the  electron  spin-destruction  rate  and  (jfE)'  1  is 
the  contribution  to  relaxation  from  spin-exchange  collisions. 
For  low  polarizations  and  small  magnetic  fields,  relaxation 
due  to  spin  exchange  is  quadratic  in  the  magnetic  field  [7] 


1  ftg  g(0)2-(27+  l)2 

W  "r^  2 


where  ^l0=BgsfjLB/ q{ 0). 

In  some  of  the  measurements  described  below,  small,  qua¬ 
sistatic  magnetic  fields  are  applied,  and  the  conditions  are 
such  that  relaxation  due  to  spin-exchange  collisions  can  be 
ignored.  In  our  present  experimental  setup  (see  Fig.  1)  opti¬ 
cal  pumping  is  along  the  z  axis  (the  longitudinal  direction) 
and  optical  rotation  of  the  probe  is  due  to  Px.  The  steady 
state  solutions  to  Eq.  (2)  can  be  found  by  setting  the  left- 
hand  side  to  zero,  and,  if  diffusion  is  neglected,  we  find 


where 


p  p  x  z  y 

x  0  b2  +  ab2  ’ 

(6) 

B2  +  A  B2 

A;' 

(7) 

AB  =  {R  +  rpl.+  I  SD )lgsP-Bi 

(8) 

Pq= sri{r  +  rpr + rSD) . 

(9) 

Note  that  Eqs.  (6)-(9)  are  independent  of  the  nuclear  spin. 

To  study  the  effects  of  spin  exchange,  we  find  it  conve¬ 
nient  to  apply  a  small  rotating  field  xBl  sin  cot+yBl  cos  wt 
in  the  presence  of  a  larger  bias  field  B„  In  this  case  we  find 
the  in-phase  and  quadrature  components  of  Px, 


p( in)  =  _  p  SsBbB  1  Am 

0  q(P )  (a>  -  no)2  +  Aar  ’ 


(10) 


n(out)  _  _  p  SsP-bB\  to  n0 

0  q(P)  (w  -  fi0)2  +  Aw2  ’ 

where  Aw=(7?  +  rpr+rSD+  II  TfE)l q(P). 

In  the  presence  of  rapid  quenching  of  the  excited  state  by 
N2,  the  effects  of  optical  pumping  and  the  optical  properties 
of  the  medium  can  be  treated  with  an  effective  ground  state 
formalism  [16].  When  the  pressure  broadened  optical  width 
is  much  larger  than  the  hyperfine  splitting,  the  optical  pump¬ 
ing  rate  for  light  of  frequency  v  is  given  by  [17] 


R  =  <3>cr  =  <S>recf 


An/2 

(An/2)2  +  (n-  n0)2  ’ 


(12) 


where  cr  is  the  absorption  cross  section,  r,,=2.8X  10-13  cm 
is  the  classical  radius  of  the  electron,  c  is  the  speed  of  light, 
/is  the  oscillator  strength  (roughly  1/3  for  7)1  light  and  2/3 
for  7)2  light),  d>  is  the  photon  flux  per  unit  area  and  An  is  the 
full  width  at  half  maximum  of  the  optical  transition  of  fre¬ 
quency  n0.  Equation  (12)  assumes  light  with  bandwidth 
much  less  than  the  pressure  broadened  optical  transition, 
valid  for  our  experimental  conditions.  Optical  rotation  of  lin¬ 
early  polarized  D2  light,  propagating  in  the  x  direction  is 
dispersive  in  the  detuning  of  the  probe  beam  from  optical 
resonance  [16,18], 


4 '>=\lrecfnP3P(v ),  (13) 

where  /  is  the  optical  path  length  and  7)(n)  =  (n-n0)/ 
[(n-n0)2  +  (An/2)2]. 


III.  EXPERIMENTAL  SETUP  AND  PROCEDURE 

The  experimental  setup  is  shown  in  Fig.  1.  A  glass  cell 
containing  a  droplet  of  Cs  metal,  600  Torr  He  buffer  gas  (to 
reduce  the  rate  at  which  atoms  in  the  central  part  of  the  cell 
diffuse  to  the  cell  walls)  and  20  Torr  of  N2  (to  eliminate 
radiation  trapping  and  improve  optical  pumping  efficiency) 
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is  placed  inside  a  four-layer  set  of  magnetic  shields.  The  cell 
has  a  roughly  cubic  profile,  about  2  cm  on  a  side.  From  the 
known  rates  of  pressure  broadening  of  Cs  lines  by  helium 
[19],  we  extrapolate  the  FWHM  of  the  D\  and  1)2  optical 
resonances  to  be  A  v=  1 5.7  and  14.1  GHz,  respectively.  The 
cell  was  heated  to  103  °C  (where  the  saturated  Cs  vapor 
concentration  is  about  [Cs]=  1.7  X  1013  cm-3)  by  flowing  hot 
air  through  the  space  between  the  walls  of  a  double-wall 
oven.  The  oven  was  designed  so  that  the  optical  path  was 
unperturbed  by  the  flowing  air. 

The  inner  magnetic  shield  has  a  cubic  shape  with  edges  of 
length  46  cm.  Magnetic  fields  were  generated  by  a  system  of 
3  coils,  wound  around  a  cubic  box  with  symmetry  axes  or¬ 
thogonal  to  each  other.  Including  image  currents  results  in 
“infinite”  solenoids  with  square  cross  sections  in  three  or¬ 
thogonal  directions.  After  degaussing,  the  residual  fields  in¬ 
side  the  magnetic  shields  are  on  the  order  of  2-3  /xG. 

Optical  pumping  was  accomplished  by  circularly  polar¬ 
ized  laser  light  propagating  in  the  z  direction  tuned  to  the 
center  of  the  Cs  D\  line  (the  exact  tuning  was  chosen  to 
minimize  light  shifts).  The  pump  beam  was  about  4  mm  in 
diameter  and  originated  from  a  Sacher  Lion  external  cavity 
diode  laser  (ECDL)  system.  The  linearly  polarized  probe 
beam,  propagating  in  the  x  direction,  had  a  cross  section 
—2X3  mm2,  and  was  tuned  about  5  optical  line  widths  from 
the  center  of  the  pressure  broadened  D2  line  (where  the  sig¬ 
nal  was  maximized).  Probe  light  was  generated  by  a  New¬ 
port  2010  laser  system.  Note  that  based  on  Eq.  (13),  one 
expects  the  maximum  optical  rotation  to  occur  for  detuning 
from  resonance  by  An/ 2.  However,  the  signal  is  the  product 
of  the  transmission  and  optical  rotation,  and  since  our  cell  is 
optically  thick  (on  resonance,  2  cm  corresponds  to  24  expo¬ 
nential  attenuation  lengths  or  an  optical  depth  (OD)0=24  un¬ 
der  the  conditions  of  our  measurements),  it  was  necessary  to 
detune  the  laser  far  from  resonance.  The  linewidth  of  both 
pump  and  probe  lasers  was  on  the  order  of  1  MHz. 

Circular  birefringence  of  the  medium  proportional  to  Px 
rotates  the  polarization  of  the  probe  beam,  which  is  analyzed 
after  the  cell  with  a  balanced  polarimeter.  To  investigate  the 
zero-field  resonance  described  by  Eq.  (6),  optical  rotation  of 
the  probe  beam  was  measured  as  a  function  of  a  static  field 
Bv  for  all  other  fields  zeroed.  Bx  (Bz)  can  be  zeroed  by  mak¬ 
ing  use  of  Eq.  (6):  a  small,  slowly  oscillating  field  (Bx)  B.  is 
applied  while  Bx  (B  )  is  adjusted  until  the  resulting  signal  is 
zero.  To  study  the  effects  of  spin-exchange  broadening  at 
nonzero  fields,  optical  rotation  was  detected  synchronously 
using  a  lock-in  amplifier  while  a  small  rotating  magnetic 
field  y Bx  cos  cot+iLBl  sin  cot  was  applied  in  the  presence  of  a 
larger  bias  field  Bz. 

IV.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 
A.  Zero-field  resonance 

In  Fig.  2  the  solid  curves  show  the  magneto-optical  rota¬ 
tion  of  a  weak  probe  beam  (7pr  —  0.5  mW/cm2)  as  a  function 
of  magnetic  field  Bx  for  two  different  pump  intensities.  The 
data  for  weak  pump  light  /pump=0.8  mW/cm2  is  well  de¬ 
scribed  by  Eq.  (6)  with  AT?  =  23  fiG.  For  more  intense  pump 
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FIG.  2.  (Color  online)  Optical  rotation  (solid  lines)  of  a  weak 
(—0.5  mW/cm2)  probe  beam  as  a  function  of  Bv  for  pump  inten¬ 
sities  as  indicated  next  to  each  trace.  The  dashed  line  overlaying  the 
high  pump  intensity  data  is  a  fit  to  two  dispersive  Lorentzians,  and 
the  dotted  line  is  a  fit  to  a  model  that  includes  the  effects  of  diffu¬ 
sion,  as  described  in  the  text. 

light  /pump=  100  mW/cm2,  optical  rotation  is  well  described 
by  the  sum  of  two  dispersive  Lorentzians  with  widths  A B 
=  56  /jlG  and  A// =  940  /xG,  as  indicated  by  the  dashed  red 
line.  In  Fig.  3  we  plot  the  width  (a)  and  peak-to-peak  ampli¬ 
tude  (b)  of  the  optical  rotation  for  the  single  feature  observed 
for  pump  intensities  below  about  15  mW/cm2  (stars)  and  for 
the  nested  features  observed  at  higher  pump  intensities 
(squares  and  triangles).  In  either  regime,  the  width  is  linear 
in  light  intensity.  When  the  two  features  become  resolved, 
the  amplitude  of  the  broad  resonance  appears  saturated, 
while  the  amplitude  of  the  narrow  resonance  continues  to 
grow,  approaching  saturation  at  the  highest  light  power. 

The  appearance  of  two  resonances  at  high  light  power  is 
essentially  due  to  diffusion  of  atoms  into  and  out  of  the 
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Pump  intensity  (mW/cm2) 

FIG.  3.  (a)  Half  width  at  half  maximum  and  (b)  peak-to-peak 
amplitude  of  the  broad,  narrow,  and  unresolved  features  observed  in 
the  zero-field  resonance  shown  in  Fig.  2,  as  a  function  of  pump 
power.  The  lines  overlaying  the  data  in  (a)  are  fits  described  in  the 
text.  The  solid  and  dashed  lines  in  (b)  serve  to  guide  the  eye. 
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pump  beam.  Similar  narrowing  effects  due  to  diffusion  have 
been  observed  in  the  context  of  electromagnetically  induced 
transparency  in  buffer  gas  cells  [21].  Nested  resonances  have 
also  been  observed  in  the  context  of  nonlinear  magneto¬ 
optical  rotation  in  paraffin  coated  cells  [22,23].  The  presence 
of  diffusion  results  in  a  layer  of  polarized  atoms  with  thick¬ 
ness  ~  \qlMY SD,  with  polarization  and  rate  of  power  broad¬ 
ening  differing  from  those  inside  the  pump  beam.  When  R  is 
small  compared  to  T SD,  the  extra  broadening  inside  the  pump 
beam  is  negligible  and  the  contributions  to  magneto-optical 
rotation  from  atoms  inside  and  outside  the  pump  beam  are 
similar,  yielding  a  single  feature.  As  R  becomes  large  com¬ 
pared  to  T SD,  the  region  inside  the  pump  beam  contributes  to 
the  broad  part  of  the  magneto-optical  rotation  curve,  while 
the  region  outside  the  beam  suffers  much  less  broadening, 
resulting  in  the  narrow  part  of  the  magneto-optical  rotation 
curve.  A  rigorous  description  of  the  nested  resonance  can  be 
obtained  by  including  the  contribution  from  diffusion  to  Eq. 
(2).  In  general,  the  solutions  are  quite  complicated,  however, 
by  neglecting  the  polarization  dependence  of  the  slowing- 
down  factor,  analytical  solutions  can  be  found  for  the  case  of 
one-dimensional  diffusion  in  Cartesian  coordinates,  valid  in 
the  limit  that  the  probe  beam  is  considerably  smaller  than  the 
pump  beam.  Figure  2  shows  a  fit  to  the  analytical  solution, 
with  the  beam  diameter  d= 2.6  mm,  qD=  15  cm2/s,  R 
=  18  000  s-\  rSD  =  465  s  1  as  free  parameters.  From  previ¬ 
ous  measurements  of  the  diffusion  constant  [20],  we  estimate 
D  =  0.46  cm2/s  for  the  buffer  gas  content  of  our  cell,  so  qD 
falls  somewhere  between  3.4  cm2/ s  and  10.1  cm2/ s  for  q  in 
the  range  of  8  and  22.  The  pump  rate  R  and  the  spin- 
destruction  rate  TSD  are  within  a  factor  of  2  of  the  values 
obtained  from  the  low  light-power  behavior  (see  below). 
Given  that  the  model  neglects  the  polarization  dependence  of 
the  slowing-down  factor,  we  consider  this  relatively  good 
agreement.  We  suspect  that  a  full  solution  taking  into  ac¬ 
count  the  polarization  dependence  of  the  slowing-down  fac¬ 
tor  will  yield  greater  accuracy  of  these  basic  parameters. 

From  the  low  light-power  data,  we  can  obtain  an  indepen¬ 
dent  measurement  of  the  intrinsic  relaxation  rate  T SD.  Over¬ 
laying  the  low  light-power  data  in  Fig.  3(a)  is  a  linear  fit 
based  on  Eq.  (8)  with  R  =  r;/punlp  yielding  r; 
=  93  s"'/(mW/cm2)  and  zero  light-power  width  A B0 
=  17 ± 3  /U.G  corresponding  to  rSD  =  300±52  s_1.  This  can 
be  compared  to  the  expected  spin-destruction  rate  based  on 
previous  measurements  of  the  spin-destruction  cross  sections 

rSD  =  [CsJ^Vg,  +  [He]FHV^  +  [N2]FN24y  (14) 

Here,  <7gD  =  “*  1CT16  cm2  [8],  4d=3X  10-23  cm2  [20],  and 
4fi=6X  10"22  cm2  [20]  are  the  spin-destruction  cross  sec¬ 
tions  for  Cs-Cs,  Cs-He,  and  Cs-N2  collisions,  respectively. 
The  mean  relative  velocity  vx  differs  between  colliding  pairs, 
and  hence  the  superscript  in  Eq.  (14).  The  contributions  to 
the  spin-destruction  rate  from  Cs,  He,  and  N2  collisions  are 
119  s-1,  88  s-1,  and  32  s_1,  respectively,  yielding  a  total 
spin-destruction  rate  rSD  =  240  s_1,  in  reasonable  agreement 
with  the  present  measurements. 


FIG.  4.  Response  of  magnetometer  to  a  small  rotating  magnetic 
field  of  magnitude  £[  =  0.55  /j,G  transverse  to  a  larger  bias  field  Bz 
for  several  different  values  of  Bz.  For  these  data  the  pump  and  probe 
intensities  were  4  mW/cm2  and  1.3  mW/cm2,  respectively. 

B.  Spin-exchange  effects 

To  explore  the  effects  of  spin  exchange  on  transverse  re¬ 
laxation  we  apply  a  bias  magnetic  field  in  the  z  direction 
(along  the  pump  beam)  and  a  small  transverse  rotating  mag¬ 
netic  field  to  excite  a  component  of  polarization  transverse  to 
the  bias  field.  In-phase  and  quadrature  components  of  the 
resulting  optical  rotation  signal  are  detected  synchronously 
using  a  lock-in  amplifier.  In  Fig.  4  we  show  the  quadrature 
sum  of  the  in-  and  out-of-phase  optical  rotation  signals  as  a 
function  of  frequency  for  several  different  values  of  the  bias 
magnetic  field.  Overlaying  the  data  are  fits  to 
aAcol  ^((o-D,0)2  +  Aar  [see  Eqs.  (10)  and  (11)].  For  these 
data,  the  pump  and  probe  intensity  were  about  4  and 
1.3  mW/cm2,  respectively.  Based  on  the  data  shown  in  Fig. 
3,  these  intensities  produce  power  broadening  by  about  a 
factor  of  2  over  the  zero  light-power  width,  however,  the 
slowing-down  factor,  determined  from  a  linear  fit  to  H0,  was 
very  nearly  q  =  22  indicating  that  the  polarization  was  quite 
low  [see  Eq.  (3)]. 

In  Fig.  5  we  plot  the  half  width  at  half  maximum  Aa>  of 
the  resonances  shown  in  Fig.  4  as  a  function  of  the  bias  field. 
Overlaying  the  data  is  a  fit  based  on  Eq.  (5)  with  Au> 
=  [?(0)T2]-1  and  D,0=gs/jLBBI q( 0),  allowing  for  a  constant 
offset  due  to  spin-destruction  collisions,  diffusion  and  power 
broadening,  yielding  a  spin-exchange  rate  1 SE 
=  14  300±350  s-1.  For  n=1.7X1013  cm-3  obtained  from 
the  saturated  vapor  pressure  curve,  Eq.  (1)  gives  a  spin- 
exchange  rate  of  about  12  000  s_l.  Temperature  fluctuations 
of  two  or  three  degrees  could  cause  significant  variations  in 
the  vapor  pressure,  and  hence  we  consider  this  reasonable 
agreement  with  measurements  of  spin-exchange  broadening. 

C.  Sensitivity 

We  evaluate  the  performance  of  the  magnetometer  by 
monitoring  the  noise  level  at  the  output  of  the  balanced  po- 
larimeter  using  a  Stanford  Research  Systems  SR770  spec- 
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Bz  (mG) 

FIG.  5.  Half  width  at  half  maximum  of  the  bias-field  resonances 
shown  in  Fig.  4  as  a  function  of  magnetic  field. 

tram  analyzer.  To  calibrate  the  magnetometer,  we  apply  a 
small  oscillating  field  By  =  Bi  cos  cot  with  Z?!  =0.55  /rG  at 
several  different  frequencies.  The  resulting  spectra  are  shown 
as  solid  lines  in  Fig.  6.  Each  trace  in  Fig.  6  is  the  result  of 
averaging  ten  spectra  obtained  during  1  s  intervals,  yielding 
a  spectral  resolution  of  1  Hz.  The  triangular  shape  of  the 
calibration  peaks  is  due  to  the  use  of  the  built-in  Hann  win¬ 
dowing  function  with  finite  spectral  resolution.  For  these 
data,  the  pump  and  probe  intensities  were  100  mW/cm2  and 
4  mW/cm2,  respectively,  and  all  three  components  of  the  dc 
magnetic_  field  have  been  zeroed.  The  sensitivity,  in 
Grins/  vHz  is  determined  by  SB  =  Bi/( \  2. S'  /  N)  yielding  a  sen¬ 
sitivity  of  about  400  pGrms/\Hz  at  30  Hz  and  about 
600  pGrms/  VHz  at  10  Hz.  The  dashed  line  in  Fig.  6  repre¬ 
sents  the  estimated  photon  shot  noise  limit  in  the  difference 
of  the  photocurrents  for  unit  bandwidth  <5/=V4e/,  where  I  is 
the  photocurrent  in  one  channel  of  the  balanced  polarimeter, 
yielding  a  photon  shot  noise  limited  sensitivity  of  about 
40  pG/VHz  at  10  Hz.  Several  possible  sources  of  the  extra 


FIG.  6.  Fourier  transform  of  magnetometer  signal  (solid  lines) 
with  calibration  peaks  of  amplitude  6[  =  0.55  fiG  applied  at  several 
different  frequencies.  The  dashed  line  represents  photon  shot  noise. 
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noise  are  (1)  real  fluctuations  of  the  ambient  field  (either  due 
to  imperfect  magnetic  shielding  or  noise  in  the  current 
source),  (2)  fluctuations  in  the  pump  power  coupled  with 
misalignment  of  the  pump  and  probe  beams,  or  (3)  vibrations 
of  the  probe  relative  to  the  pump  in  the  plane  common  to 
both  beams.  Optimization  of  geometry  to  maximize  the  over¬ 
lapping  volumes  of  the  pump  and  probe  beams  will  likely 
yield  improvements  in  the  photon  shot  noise  limit. 

V.  CONCLUSIONS 

We  demonstrated  a  Cs  atomic  magnetometer  in  the  spin- 
exchange-relaxation-free  regime.  The  primary  advantage  of 
using  Cs  is  the  ability  to  work  at  lower  temperatures.  Future 
work  with  atomic  magnetometers  in  the  context  of  microflu¬ 
idic  NMR  will  make  use  of  this  feature.  At  103  °C  we  real¬ 
ized  magneto-optical  rotation  features  with  intrinsic  line- 
widths  of  17  piG  corresponding  to  a  relaxation  rate  of  about 
300  s-1  when  the  spin-exchange  rate  was  about  rSP 
=  14  000  s"1.  We  achieved  a  sensitivity  of  400  pG/\Hz. 
Based  on  estimates  of  the  photon  shot  noise,  we  project  a 
sensitivity  of  about  40  pG/  VHz.  We  suspect  that  the  demon¬ 
strated  sensitivity  was  limited  by  pump  laser  noise  and  am¬ 
bient  magnetic  field  noise.  Theoretical  optimization  of  the 
magnetometer  (presented  in  the  Appendix  below)  indicates  it 
should  be_  possible  to  achieve  sensitivity  on  the  order  of 
2  pG/  VHz  in  a  1  cm3  volume.  We  believe  the  primary  rea¬ 
son  for  falling  short  of  this  level  is  suboptimal  geometry 
(probe  beam  cross  section  was  only  2X3  mm2)  detuning 
and  light  power  of  both  pump  and  probe. 
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APPENDIX:  THEORETICAL  OPTIMIZATION 

We  now  present  a  theoretical  optimization  of  the  magne¬ 
tometer,  maximizing  sensitivity  to  small,  quasistatic  fields. 
The  analysis  is  similar  to  Refs.  [24,25],  in  that  spin- 
projection  noise  and  photon  shot  noise  are  considered  inde¬ 
pendently  (noise  due  to  light  shifts  is  not  considered,  because 
in  principle,  it  can  be  eliminated  by  orthogonality  of  pump 
and  probe  beams  [25]).  Spin-projection  noise  is  typically 
written  as 

F=,  (Al) 

yyNtT2 

where  y  is  the  gyromagnetic  ratio,  N  is  the  number  of  atoms, 
t  is  the  measurement  time,  and  T2  is  the  transverse  relaxation 
time.  In  the  SERF  regime  T2  =  q(P)/TSD  and  y=qxfiHl  q(P) 
both  depend  on  the  nuclear  slowing-down  factor.  Inserting 
these  expressions  into  Eq.  (Al),  one  might  conclude  that  the 
atomic  shot  noise  limit  scales  as  \jq(P).  However,  it  turns  out 
that  the  nuclear  slowing-down  factor  drops  out  of  the  prob¬ 
lem.  The  reasons  are  somewhat  subtle,  so  we  go  into  some 
detail. 

Spin-projection  noise  arises  (in  the  present  geometry)  due 
to  uncertainty  in  the  x  component  of  angular  momentum  Fx, 
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defined  as  A Fx=  \j(Fx)-(Fx)2.  For  a  spin-temperature  distri¬ 
bution  with  polarization  in  the  z  direction  poce0G  where  f3 
=  ln(l+P)/(l-P)  is  the  spin-temperature  parameter  [11] 
(Fx)  =  0.  Thus  A Fx=  \  Tdppf).  Evaluation  of  this  trace  results 
in  A/-’t(/J)=  \q(P)/4  per  atom  with  p  normalized  so  that 
Trp=l.  Assuming  that  N=nV  uncorrelated  atoms  are  in¬ 
volved  in  the  measurement,  the  ensemble  averaged  uncer¬ 
tainty  scales  as  1  /  \fn  V, 


In  the  large  polarization  limit,  <SFX(1)  =  y2/nV.  This  limit  can 
be  obtained  from  the  angular  momentum  commutation  rela¬ 
tions  \Fx,Fy\=iFv  which  yield  the  minimum  uncertainty 
VK^>|/2.  If  all  the  atoms  are  in  the  stretched  state,  corre¬ 
sponding  to  P=  1,  (F.)  =  4,  we  have,  again  assuming  uncor¬ 
related  atoms,  SFx=i2l nV.  The  uncertainty  in  the  low  po¬ 
larization  limit  is  somewhat  larger,  SFx(0)  =  ^HI2nV.  This 
limit  can  also  be  verified  by  noting  that  for  an  unpolarized 
sample  p=  1/(25+  1)(2/+ 1)  and  Tr(pF2J=Tr(pF2).  As  an 
aside,  we  note  that  the  reduction  in  uncertainty  SFX  with 
increasing  polarization  only  occurs  for  angular  momentum 
greater  than  1/2. 

After  measuring  continuously  for  time  t  long  compared  to 
the  lifetime  of  the  polarization  q(P)l  (R  +  Y  fr+Y SD),  the  un¬ 
certainty  is  [24,26] 


(SFX),=  SFX 


2  q(P) 


(R  + 1  Pi-+ r$D)f 


(A3) 


q(P) 

V 2t(R  +  rpr + r  SD)n  V 


(A4) 


In  a  spin-temperature  distribution,  the  ratio  of  the  total  angu¬ 
lar  momentum  to  that  stored  in  the  electron  is  given  by  q{P ), 
<F)  =  [p(F)/2]P  [11]  and  thus 


SPr 


q{P ) 


< SFX), 


t{R  +  rpr + rSD)nv 


(A5) 


The  uncertainty  SBy  in  a  measurement  of  By  is  related  to 
flucutations  of  Px  via  Eq.  (6)  (for  all  other  fields  zeroed), 


SB, 


R  +  I  pi  +  I  sd  8PX 
gsPB  PZ 


(A6) 


Inserting  Eq.  (A5)  into  Eq.  (A6)  we  find  that  the  spin- 
projection  noise  is 


i  J2(r  +  rpr+ rSD) 

gsPRPz.  '  n 


(A7) 


It  is  interesting  to  note  that  this  result  is  independent  of  any 
nuclear  slowing-down  factors.  Neglecting  broadening  due  to 
the  probe  beam,  the  minimum  value  of  spin-projection  noise 


SB 


min 

spn 


is  obtained  when  /\’  =  2FSD. 


We  now  address  photon  shot  noise.  To  simplify  the  analy¬ 
sis,  we  assume  that  the  volume  V  occupied  by  the  sample  is 
a  cube  with  sides  of  length  /,  fully  illuminated  by  both  pump 
and  probe.  If  the  probe  beam  is  detuned  far  from  resonance 
so  that  the  medium  is  optically  thin,  optical  rotation  drops 
slowly,  scaling  as  D{v)~  ll(v-v0),  compared  to  absorption 
which  scales  as  1/(v-Vq)2,  a  very  favorable  situation.  In  this 
case,  photon  shot  noise  in  the  optical  rotation  angle  is  given 
by  S(/)=  1  /  2  v(fi0/2t,  where  <I»0  is  the  probe  photon  intensity 
and  / 2  is  the  cross  section  of  the  probe  beam.  Combining  this 
with  Eqs.  (13)  and  (A6),  the  photon  shot  noise  contribution 
to  magnetic  field  sensitivity  is 


SBn 


i  2(r  +  rpr + 1 SD) 


psn  gspBPz  lrecfnD(v) V<3>o l2t 
This  can  be  rearranged 

2(r  + 1  p,  +  rSD) 


SB 


1 


psn  ' 


(A9) 


(A  10) 


gsPBPz^nVt  Vrpr(OD)0 

where  (OD)0=2recfnII  Av  is  the  optical  depth  on  resonance 
and 


rpr  =  $0  recf 


Avl2 
(v-  v0 )2 


(All) 


is  the  probe  rate  for  far  detuned  light. 

Adding  spin-projection  noise  Eq.  (A7)  and  photon  shot 
noise  Eq.  (A  10)  in  quadrature,  yields 


SB : 


1 


gsFuP  z^nVt 


2  (r  + 1  pr+ rSD)  + 


4  (r  +  r  pr + rSD) 


rpr(OD)0 


(A  12) 


For  a  cubic  volume  of  1  cm3  and  a  density  of  n=1.7 
X  1013  cm-3,  the  resonant  optical  depth  (OD)0~  12  for  the 
buffer  gas  contents  of  our  cell.  Assuming  a  spin-destruction 
rate  rSD=300  s-1  obtained  in  the  experiment  at  103  °C,  Eq. 
(A12)  reaches  a  minimum  of  about  2.4  pG/ \Hz  with  pump 
and  probe  rates  R  =  7 10  s~!  and  rpr=91  s~'  (here  we  have 
assumed  that  a  bandwidth  of  1  Hz  corresponds  to  a  measure¬ 
ment  time  of  0.5  s). 

Further  inspection  of  Eq.  (12)  shows  that  the  second  term 
underneath  the  radical  can  be  made  small  if  the  resonant 
optical  depth  is  large.  In  this  case,  rpr  can  be  made  small  so 
that  the  first  term  underneath  the  radical  is  minimized.  In  the 
limit  of  infinite  resonant  optical  depth  Eq.  (A  12)  is  optimized 
for  R  =  2Tsd  and  rpr~0,  in  which  case  Eq.  (A  12)  reduces  to 
Eq.  (A8).  For  a  volume  of  1  cm3  and  a  density  of  »=1.7 
X  1013  cm-3,  and  spin-destruction  rate  rSD=300  s~*  ob¬ 
tained  in  the  experiment  at  103  °C,  Eq.  (A8)  evaluates  to 
about  1.2  pG/  VHz.  Slightly  better  sensitivity  may  be 
achieved  by  operating  at  higher  densities  where  Cs-Cs  colli¬ 
sions  completely  dominate  the  collisional  spin-destruction 
rate  [see  the  discussion  following  Eq.  (14)]. 

The  above  analysis  indicates  that  optimal  sensitivity  is 
achieved  when  the  probe  is  tuned  sufficiently  far  from  reso¬ 
nance  so  that  the  medium  is  optically  thin,  minimizing  pho¬ 
ton  shot  noise.  However,  technical  sources  of  noise  due  to, 
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e.g.,  vibrations  or  air  currents,  are  often  much  larger  than 
photon  shot  noise.  In  this  case,  it  is  desirable  to  tune  the  laser 
closer  to  optical  resonance  so  that  the  optical  rotation  due  to 
small  magnetic  fields  is  larger  than  other  sources  of  noise. 
Finally,  we  note  that  the  optimal  tuning  of  the  probe  light 
depends  on  the  particulars  of  the  probing  scheme.  For  ex¬ 
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ample.  Ref.  [27]  considers  the  case  of  nonlinear  magneto¬ 
optical  rotation,  where  optimal  sensitivity  is  achieved  when 
the  probe  is  tuned  so  that  there  is  roughly  one  optical  depth. 
The  primary  reason  for  the  difference  is  that  optical  rotation 
in  that  case  scales  (similarly  to  absorption)  as  the  inverse 
square  of  the  detuning. 
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